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SUMMARY 


The  spanwise  lift  distributions  of  -wings  of  low  aspect  ratio  but  of 
arbitrary  plan  form  and  angle -of -attack  distribution  have  been  analyzed 
by  two  well-established  concepts:  the  virtual -mass  concept  and  the 
Weissinger  method <>  Both  concepts  have  been  found  to  yield  the  same 
simple  integral  expression  for  the  spanwise  lift  distribution  in  terms 
of  the  spanwise  angle -of -attack  distribution.  Tables  and  figures  of 
lift  distributions  of  low-aspect -ratio  wings  have  been  presented. 

Within  certain  limitations  these  distributions  are  independent  of  the 
plan  form.  For  plan  forms  with  leading  edges  swept  back  behind  the  Mach 
cone,  the  results  of  the  analysis  are  applicable  at  supersonic  speeds. 


INTRODUCTION 


For  wings  of  very  low  aspect  ratio,  of  the  order  of  1  or  less,  at 
small  angles  of  attack  the  flow  may  be  assumed  to  be  esentially  two- 
dimensional  in  planes  normal  to  the  stream  direction.  (See  reference  1.) 
In  the  present  paper  this  concept  is'  shown  to  lead  to  a  simple  formula¬ 
tion  of  the  problem  of  calculating  the  spanwise  lift  distribution  of 
wings  of  very  low  aspect  ratio  for  various  angle -of -attack  distributions 
in  terms  of  conditions  at  the  widest  section  of  the  wing.  The  same 
problem  is  also  formulated  as  a  limiting  case  of  the  Weissinger 
L-method*  (See  reference  2.)  Some  comparisons  are  drawn  between  the 
assumptions  and  limitations  of  the  two  approaches. 

Spanwise  lift  distributions  have  been  calculated  for  angle -of - 
attack  distributions  which  correspond  to  several  symmetrical  and 
antisymmetrical  twists,  as  well  as  to  the  deflection  of  flaps  and 
ailerons  of  various  spans.  The  results  are  presented  in  a  table  and  in 
several  figures. 


SYMBOLS 


A 

b 

>a 

bf 

c 


c* 


c 


aspect  ratio 

wing  span  (at  widest  section) 
aileron  span  (total  for  both  wings) 


flap  span  (total  for  both  wings) 
chord 


dimensionless  chord 


average  chord 


cl 

CL 

CLl/2 

CDi 

C7 

CBM 

L 

Ll/2 

M 

q. 

s 

v 


local  lift  coefficient 


wing  lift  coefficient 


lift  coefficient  on  either 


half  of  wing 


induced  drag  coefficient 

rolling-moment  coefficient  (Rollingoment^ 

bejading -moment  coefficient 

lift 

lift  on  either  half  of  wing 


vertical  momentum  of  virtual  mass  per  unit  length 
dynamic  pressure 
wing  area 
airspeed 

downwash  velocity 
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y 

y* 

y* 

^o* 

a 

a6 

r 

r* 

7 

6 

Tl* 

Q 

<6 

A 


spanwis©  distance 

dimensionless  spanwise  distance  (fraction  of  semispan) 

dimensionless  spanwise  distance  of  center  of  pressure  (fraction 
of  semispan)  / 

dimensionless  spanvise  distance  at  point  of  discontinuity  in 
angle  of  attack 

angle  of  attack  (slope  of  mean-camber  surface) ,  radians 
angle  of  attack  at  tip,  radians 


flap  or  aileron  effectiveness  parameter 


?ter  ^ 


boi/bb 

be  1  /ba 


circulation 


dimensionless  circulation  (—  or  c*ct  ) 

VbV 

loading  coefficient  or 


flap  or  aileron  deflection,  radians 
variable  of  integration  corresponding  to  y* 
trigonometric  variable  corresponding  to  y*(cos“^  y*) 
variable  of  integration  corresponding  to  Q 
sweep  angle 


ANALYSIS 

The  Virtual-Mass  Approach 


According  to  the  theory  of  reference  1  the  disturbance  set  up  by 
wings  of  low  aspect  ratio  at  small  angles  of  attack  is  essentially  two- 
dimensional  in  planes  perpendicular  to  the  free  stream.  The  lift  on  a 
spanwise  strip  of  width  dx  may  then  be  obtained  from  the  time  rate  of 


k 
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change  of  the  vertical  momentum,  of  the  virtual  mass  associated  with  the 
given  two-dimensional  flow 


dL  =  M  doc 
dt 

=  v™  dx 
dx 


=  Y  dM 

where  M  is  the  vertical  momentum  of  the  virtual  mass  per  unit  length. 
The"  lift  on  the  entire  wing  is  then  . 

L  =  VMW  (D 

■where  is  the  downward  vertical  momentum  at  the  maximum  width  of 

the  airfoil.  Any  portion  of  the  wing  "behind  the  section  of  maximum 
width  is  assumed  to  "be  alined  with  the  wake  so  that  it  does  not  affect 
the  lif t . 

Equation  (l)  indicates  that  the  lift  depends  entirely  on  the 
geometry  of  the  widest  section.  A  physical  interpretation  of  this 
equation  is  that  the  wake  leaves  the  airfoil  parallel  to  the  surface  at 
the  section  of  maximum  width  and  that  the  vorticity  distribution  in  the 
wake  is  such  as  to  produce  the  required  downwash  distribution  in  any 
plane  perpendicular  to  the  wake.  The  determination  of  the  spanwise  lift 
distribution  on  the  wing  then  reduces  to  a  problem  almost  identical  in 
form  with  that  of  two-dimensional  thin-airfoil  theory  and  very  similar 
to  but  simpler  than  the  classical  lifting-line  theory. 

If  the  bound  circulation  on  the  wing  at  any  spanwise  location  is 
given  by  the  parameter 


7  = 


and  if  the  spanwise  distance  is  given  by  the  parameter 
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and  an  associated  variable  of  integration  tj*,  then  the  down-wash  angle 
at  any  spanwise  location  y*  in  the  wake  is  given  by 


v  _1_  f 1  d7_  dn* 

V  =  2jtA  J  dt]* ,  y*  -  T)* 

From  the  foregoing  considerations  this  angle  must  be  equal  to  a,  the 
slope  of  the  mean-camber  surface  at  the  widest  section.  If,  furthermore, 
the  trigonometrical  variables  8  and  d  are  defined  as 

9  =  cos“-L  y* 

■8  =  COS~l  Tj* 

equation  (2)  becomes 


T*  3Z.  dfl 

2«A  Jq  dd  cos  -5  -  cos  9 

Equation  (3)  may  be  solved  by  means  of  an  inversion  integral: 


7  =  =A  a  sin  ^  log 

*  J  0 


•8  +  9 


1$  -  01 


For  a  constant  angle  of  attack,  for  instance,  equation  (1)  yields  an 
elliptic  lift  distribution  with  a  lift  coefficient  equal  to  ^  Aa  in 
agreement  with  the  result  of  reference  1. 


Three -Quarter-Chord  Approach 

Another  approach  to  the  problem  of  calculating  spanwise  lift  dis¬ 
tributions  on  wings  of  low  aspect  ratio  is  to  consider  the  limit,  as  the 
aspect  ratio  goes  to  zero,  of  the  equations  of  the  Weissinger  L-method. 
This  method  consists  of  concentrating  the  bound  vorticity  at  the  quarter- 
chord  line,  calculating  the  downwash  angle  at  the  three -quarter -chord 
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line,  and  setting  that  dovnvash  angle  equal  to  the  slope  of  the  mean- 
camber  surface  at  the  three-quarter-chord  line.  This  equation  may  be 
•written  as 


_i_  r1  dT*  at)* 

dry*  JT*  -  TJ* 


ar* 

dry* 


dt)*  =  a 


vhere  the  dimensionless  circulation  r*  is  defined  as 


(5) 


r*  =  — 

1  tv 


and  the  plan-form  function  F  is  defined,  for  tj*  ^  0,  as 


F  = 

y*  _ 


and,  for  tj*  =  0,  as 


=  c*/2  J  N 


1  + 


tan  A 


2  *  ( W)2 


y*  ■  ry 


1  +  2$72  tan  A 


1 


J 


+  2  tan  A 


C1  *  tan  A)2  *  ( 


1  +  2-£j2  tan  A 


and  -where  a  now  represents  the  slope  of  the  mean -camber  surface  at 
the  three-quarter-chord  line. 

As  the  aspect  ratio  approaches  zero,  the  F  function  approaches 
the  value  tan  A,  which  is  constant  if  the  semispan  quarter-chord  line 
is  straight.  For  instance,  all  plan  forms  of  figure  1  except  that  of 
figure  l(c)  have  straight  quarter-chord  lines.  If  F  is  constant  the 
second  integral  in  equation  (5)  vanishes.  Furthermore,  the  dimen¬ 
sionless  chord  c*  approaches  infinity  and  its  reciprocal  vanishes. 
Consequently,  the  second  term  of  equation  (5)  ma y  be  expected  to 
approach  zero  quite  rapidly  as  the  aspect  ratio  approaches  zero. 
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If  the  parameter  7,  -which  is  related  to  the  dimensionless 
circulation  r*  by 

7  =  £r* 

2 

and  the  trigonometric  variables  are  substituted  in  equation  (5)  and  if 
the  second  term  is  neglected,  equation  (5)  becomes  identical  with 
equation  (3),  except  that  a  now  stands  for  the  angle  of  attack  at  the 
three-quarter-chord  line.  The  solution  for  the  circulation  is  then 
given  by  equation  (k)  as  before. 


RESDUTS 


In  order  to  facilitate  the  application  of  equation  (k)  to  cases  of 
practical  interest,  the  integral  of  that  equation  has  been  evaluated 
for  several  angle -of -attack  distributions  which  can  be  expressed 
analytically  as  follows: 

(a)  The  constant  angle -of -attack  or  additional-loading  case 
analyzed  in  reference  1 

(b)  The  linear  antisymmetrical-twist  or  damping-in-roll  case 
previously  analyzed  in  reference  3 

(c)  The  symmetrical-  and  antisymmetrical-twist-distribution 

cases  a  =  Oj.Jy*|n  in  which  ag  is  the  angle  of  attack  at  the  tip  and 
is  considered  to  have  equal  and  opposite  values  at  the  two  tips  in  the 
antisymmetrical  case 

(d)  The  flap  case  in  which  a  =  058  from  y*  =  -yD*  to  y*  =  yD* 
and  a  =  0  elsewhere,  agS  being  the  effective  angle  of  attack  due  to 
flap  deflection  (some  considerations  involved  in  the  selection  of  the 
value  of  05  are  outlined  in  the  section  entitled  '‘Discussion ") 

(e)  The  aileron  case  in  which  a  =  ±ag8  outboard  of  y*  =  ±y0* 
and  a  =  0  elsewhere 

The  resulting  lift  distributions  for  cases  (a)  to  (e)  are  presented 
in  table  1  and  in  figures  2  to  5*  All  results  are  presented  for  unit 
aspect  ratio  A  and  unit  angle  of  attack  a,  where  a  is  the  angle  of 
attack  at  the  tip  for  the  twist  cases  and  is  the  effective  angle  of 
attack  due  to  control  deflection  a§8  for  the  flap  and  aileron  cases. 
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The  twist  distributions  of  case  (c)  may  be  useful  in  the  analysis 
of  structural  twists •  If  a  structural  twist  distribution  is  of  the 
type  considered  in  case  (c) ,  the  corresponding  lift  distribution  may  be 
taken  from  the  table  or  the  figures  directly.  Otherwise,  a  structural 
twist  distribution  can  usually  be  approximated  by  the  linear  super¬ 
position  of  two  or  more  of  the  twist  distributions  covered  in  case  (c) . 
The  lift  distribution  is  then  given  by  a  superposition  of  the  lift  dis¬ 
tributions  corresponding  to  those . particular  twist  distributions • 

The  results  for  case  (d)  may  be  applied  to  the  case  of  symmetrically 
deflected  ailerons  by  subtracting  a  flap  distribution  from  the  constant- 
angle  -of  -at  tack  distribution  for  the  same  angle.  Similarly,  the  results 
for  case  (e)  may  be  used  to  calculate  the  lift  distribution  due  to 
inboard  aileron’s  by  subtracting  the  aileron  distribution  for  the  given 
value  of  y0*  from  the  aileron  distribution  for  100 -percent -semispan 
aileron. 

From  the  lift  distributions  the  lift,  rolling,  and  bending  moments, 
as  well  as  the  lateral  center  of  pressure  of  the  wing,  may  be  calculated 
in  the  form  of  dimensionless  coefficients  % 

r  ** 

Cp,  =  i  /  7  sin  6  <30 

2  J  o 


cLi/2  = 


sin  9  <3 6 


C  i  =  y  sin  26  dJ9 


i  r*/2 

CEM  =  5  J0  7  sin  26  d0 


y*  = 


CEM 

°Ll/2 
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The  values  of  these  coefficients  for  the  various  angle -of -attack 
distributions  are  given  in  table  lj  the  lift  coefficient  due  to  flap 
deflection  and  the  rolling  moment  due  to  aileron  deflection  are  also 
given  in  figure  6  as  functions  of  the  flap  or  aileron  span  ratio. 

The  induced  drag  coefficient  may  be  calculated  from  the  conventional 
lifting-line  equation  since  the  drag  is  not  affected  by  the  longitudinal 
distribution  of  the  lift.  Consequently, 

=  i  J  ay  sin  8  68  (6) 


For  the  constant  angle-of -attack  case  equation  (6)  becomes 


as  also  given  in  reference  1. 


DISCUSSION 


The  two  approaches  to  the  problem  of  calculating  the  spanwise  lift 
distributions  of  wings  of  low  aspect  ratio  presented  in  this  paper  have 
in  common  all  the  limitations  of  linear  potential  theory  and  are  subject 
to  certain  additional  inherent  limitations  as  well.  These  additional 
limitations  and  the  effect  they  have  on  the  validity  of  the  results  of 
the  two  approaches  are  discussed  in  the  succeeding  paragraphs . 

The  validity  of  the  virtual-mass  approach  has  been  demonstrated 
for  the  constant -angle -of -attack  case  in  references  1  and  4  and  for 
the  linear  antisymmetric  angle -of -attack  case  in  reference  5  primarily 
by  comparison  with  more  rigorous  analytical  methods.  The  validity  of 
the  Weissinger  L-method  for  plan  forms  of  low  aspect  ratios  has 
similarly  been  demonstrated  for  the  constant -angle -of -attack  case  in 
references  2  and  4  and  for  the  linear  antisymmetric  case  in  reference  5 
by  comparison  with  both  analytical  and  experimental  results.  The 
extension  of  the  approach  of  reference  1  and  the  specialization  of  the 
method  of  reference  2  presented  in  this  paper  may  therefore  also  be 
expected  to  be  valid  provided  certain  conditions  of  aspect  ratio,  plan 
form,  angle  of  attack,  and  Mach  number  are  met. 

The  aspect  ratio  has  been  assumed  to  be  vanishingly  small  in  the 
analysis .  For  the  purpose  of  comparison  the  lift-curve  slopes  and 
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additional  lift  distributions  calculated  by  the  method  of  reference  2 
are  shorn  in  figure  7  for  various  plan  forms  of  aspect  ratio  1.5 
together  with  the  values  calculated  ty  the  low-aspect-ratio  analysis  of 
the  present  paper,  namely  an  elliptic  lift  distribution  and  a  lift- 

curve  slope  of  ^-(1.5)  =  2.356.  At  an  aspect  ratio  of  1*5  the  lift  dis¬ 
tributions  for  the  various  plan  forms  are  in  fair  agreement  with  the 
elliptic  distribution  and  with  each  other.  The  lift-curve  slopes  of 
the  various  plan  forms  are  almost  the  same  and  are  about  15  percent 
lower  on  the  average  than  that  calculated  by  the  simplified  analysis. 

This  type  of  agreement  does  not  appear  to  be  sufficient  for  most 
practical  purposes.  In  reference  1,  however,  the  lift -curve  slope  is 
shown  to  be  overestimated  by  only  about  4  percent  at  an  aspect  ratio 
of  1  as  compared  with  a  more  exact  lifting-surface  method.  Consequently, 
an  aspect  ratio  of  1  may  in  general  be  regarded  as  the  upper  limit  at 
-which  the  results  of  the  analysis  of  this  paper  may  be  expected  to 
apply. 

The  virtual -mass  concept  utilized  in  reference  1  does  not  imply  any 
plan-form  or  angle -of -attack  restrictions  except  that  the  aspect  ratio 
must  be  low  and  that  there  cannot  be  any  abrupt  longitudinal  Changes  in 
the  span  or  the  slope  of  the  surface •  As  pointed  out  in  reference  1, 
however,  it  is  convenient  to  apply  this  concept  only  to  pointed  plan 
forms  such  as  those  shown  in  figures  l(b),  l(d),  l(e),  l(f),  l(h),  l(i), 
and  l(  j ) .  Wings  with  re-entrant  trailing  edges  cannot  be  conveniently 
treated  by  this  method.  Furthermore,  any  part  of  the  wing  downstream 
of  the  widest  section  must  be  alined  with  the  widest  section.  Any  flaps 
and  ailerons  (such  as  those  shown  in  fig.  1(h))  must  themselves  be  of  low 
aspect  ratio  in  order  for  the  low-aspect  -ratio  analysis  to  apply  to  the 
wing. 

The  Weissinger  L-method,  on  the  other  hand,  applies  to  all  plan 
forms  except  those  with  curved  or  discontinuous  quarter-chord  lineB  or 
discontinuous  chord  distributions  (it  does  not,  for  instance,  apply  to 
an  elliptic  plan  form) .  It  has  the  disadvantage,  however,  of  being 
strictly  applicable  only  to  uncambered  wings.  This  observation  results 
from  the  following  reasoning:  The  Weissinger  method  is  based  on  the 
slope  of  the  mean-camber  surface  at  the  three-quarter-chord  line  at  all 
aspect  ratios.  For  wings  of  very  low  aspect  ratio  the  virtual-mass 
concept  indicates  that  the  widest  section  determines  the  lift  on  the 
wing.  Consequently,  the  application  of  the  Weissinger  method  to  wings 
of  low  aspect  ratio  can  give  correct  results  only  if  the  wing  is 
uncambered,  so  that  the  slope  at  the  three -quarter-chord  line  is  equal 
to  that  at  the  widest  section. 

Consequently,  the  virtual -mass  concept  is  (for  practical  purposes) 
more  restricted  as  to  plan  form  but  less  restricted  as  to  camber  than 
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the  three-quarter-chord  method .  In  view  of  the  agreement  of  the  two 

*  concepts  in  the  cases  to  which  both  apply,  however,  slight  deviations 
from  the  plan-form  conditions  of  the  one  or  the  camber  conditions  of  the 
other  may  be  expected  to  have  little  effect  on  the  results. 

If  either  concept  is  used  to  analyze  a  wing  with  deflected  flaps 
or  ailerons,  the  question  arises  as  to  what  value  to  choose  for  the 
parameter  0,5.  According  to  the  virtual-mass  concept  a  particular  flap 
or  aileron  angle  will  be  just  as  effective  in  producing  lift  as  the  same 
angle  imposed  on  the  entire  section;  therefore  the  parameter  05  is  1 
and  the  flap  or  aileron  chord  ratio  does  not  enter  into  the  problem. 

On  the  other  hand,  the  three-quarter-chord  concept  Is  applicable  only  if 
the  flap  or  aileron  occupies  100  percent  of  the  chord,  in  which  case  ag 
is  again  1.  In  the  more  conventional  case  of  partial-chord  flaps  and 
ailerons  it  might  be  expected,  by  analogy  to  the  trailing-edge  concept, 
that  a  value  of  05  between  1  and  the  two-dimensional  value  of 
about  l/2  should  be  chosen  rather  than  the  two-dimensional  value  for 
the  given  flap  or  aileron  chord  ratio;  the  determination  of  the  exact 
value  would  require  either  experimentation  or  a  more  refined  theory. 

The  angle  of  attack  has  been  assumed  small  in  the  analysis  for 
several  reasons:  At  large  angles  of  attack  the  flows  in  planes  perpen¬ 
dicular  to  the  direction  of  flight  cannot  be  considered  to  be  inde¬ 
pendent  of  each  other,  the  trailing  vorticity  is  not  shed  substantially 
in  the  plane  of  the  wing  but  more  nearly  parallel  to  the  direction  of 
flight,  and,  finally,  at  large  angles  of  attack  the  flow  separates  and 
potential -flow  considerations  no  longer  apply. 

*  . 

If  the  leading  edge  and  the  sides  of  the  wing  (in  the  case  of  plan 

forms  of  the  type  shown  in  figs.  l(a),  l(b),  l(e),  and  l(f),  for 
*  instance)  are  sharp  the  flow  may  detach  at  very  low  angles  of  attack. 

The  flow  separation  affects  both  the  total  lift,  which  assumes  a 
parabolic  variation  with  angle  of  attack  rather  than  the  usual  linear 
one,  and  the  lift  distribution.  (See  reference  6.)  The  theoretical 
treatment  of  this  phenomenon  transcends  the  scope  of  this  paper;  an 
analysis  of  a  simple  case  is  presented  in  reference  7  •  If  the  leading 
edge  and  sides  are  well  rounded,  the  flow  separation  and  the  resulting 
parabolic  variation  of  the  lift  coefficient  and  distortion  of  the  lift 
distribution  tend  to  be  postponed  to  higher  angles  of  attack.  At  lower 
angles  of  attack  the  analysis  of  this  paper  may  be  applied  to  the 
problem. 

The  virtual -mass  approach  has  been  shown  to  be  valid  for  supersonic 
speeds  (reference  l),  provided  the  leading  edge  is  swept  back  suffi¬ 
ciently  far  behind  the  Mach  cone,  as  well  as  for  subsonic  speeds. 
Computation  of  the  downwash  in  the  manner  utilized  in  this  paper  is  also 
f  valid  at  supersonic  speeds  because  far  downstream  of  the  cross  element 
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of  a  horseshoe  vortex  the  down-wash  is  the  same  whether  the  vortex  is  in 
subsonic  or  supersonic  flow*  Consequently,  the  three -quarter -chord 
concept  may  also  be  expected  to  be  valid  at  supersonic  speeds  provided 
the  leading  edge  is  swept  back  far  behind  the  Mach  cone  and  provided 
the  angle -of -attack,  camber,  and  plan-form  conditions  are  such  that  the 
analysis  of  this  paper  applies. 


CONCLUDING  REMARKS 


The  spanvise  lift  distributions  of  wings  of  low  aspect  ratio  but 
of  arbitrary  plan  form  and  angle -of -attack  distribution  have  been 
analyzed  by  two  well-established  concepts;  the  virtual-mass  concept 
and  the  Weissinger  method.  Both  concepts  have  been  found  to  yield  the 
same  simple  integral  expression  for  the  spanwise  lift  distribution  in 
terms  of  the  spanwise  angle-of -attack  distribution.  Tables  and  figures 
of  lift  distributions  of  low-aspect-ratio  wings  have  been  presented. 
Within  certain  limitations  these  distributions  are  independent  of  the 
plan  form.  For  plan  forms  with  leading  edges  swept  back  far  behind 
the  Mach  cone,  the  results  of  the  analysis  are  applicable  at  supersonic 
speeds . 


Langley  Aeronautical  Laboratory 
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Langley  Air  Force  Base,  Va.,  August  25,  19^9 
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Figure,  Dr  Fft  o//sfnbut/on .  of  tv  mgs  of  low  aspect  ratio  for  constant 
angle  of  attacK  and  for  symmetrical  twists. 
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Figure  3 —  Lift  distributions  of  wings  of  low  aspect  ratio  for  antisymmetries!  twists. 


Figure  5*—  Lift  distributions  of  vings 


Rot l/np-moment  coefficient  Lift  coefficient  due 

due  to  a/teron  deflection >,  fl  to  f/op  deflection, 
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Figure  6.—  Lift  coefficient  due  to  flap  defxection  and  rolling— moment 
coefficients  due  to  aileron  deflections  for  wings  of  low  aspect 
ratio. 


